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Abstract. Over the past decade measurements of unpolarized structure functions with 
unprecedented precision have significantly advanced our knowledge of nucleon structure. These 
have for the first time allowed quantitative tests of the phenomenon of quark-hadron duality, 
and provided a deeper understanding of the transition from hadron to quark degrees of 
freedom in inclusive scattering. Dedicated Rosenbluth-separation experiments have yielded 
high-precision transverse and longitudinal structure functions in regions previously unexplored, 
and new techniques have enabled the first glimpses of the structure of the free neutron, without 
contamination from nuclear effects. 



1. Introduction 

Throughout the modern history of nuclear and particle physics, measurements of structure 
functions in high-energy lepton-nucleon scattering have played a pivotal role. The demonstration 
of structure function scaling in early deep-inelastic scattering (DIS) experiments at SLAC in 
the late 1960's [1] established the reality of quarks as elementary constituents of protons and 
neutrons — a feat recognized by the award of the 1990 Nobel Prize in Physics to Friedman, 
Kendall and Taylor. This paved the way to the development of quantum chromodynamics 
(QCD) as the theory of strong nuclear interactions, and its subsequent confirmation several 
years later through the discovery of logarithmic scaling violations in structure functions [2] . 

The interpretation of structure functions in terms of quark and gluon (or parton) momentum 
distributions resulted in the emergence of a remarkably simple and intuitive picture of the 
nucleon [3], allowing a vast amount of scattering data to be described in terms of a few universal 
functions — the parton distribution functions (PDFs). At leading order (LO) in a^, the F2 
structure function of the proton, for example, could be simply represented as a charge squared- 
weighted sum of PDFs, 

F2 = Y.4 <^ + 9) = + ^) + \<d + d) + ^x{s + s) + • • • , (1) 

where q and q are the quark and antiquark momentum distribution functions, usually expressed 
as functions of the momentum fraction x of the nucleon carried by the parton, at a scale given 
by the momentum transfer squared Q^. 

Over the ensuing decades concerted experimental DIS programs at SLAC, CERN, DESY 
and Fermilab have provided a detailed mapping of the PDFs over a large range of kinematics, 
with and x spanning several orders of magnitude. To manage the ever increasing number 
of data sets, from not only inclusive DIS but also other high energy processes such as Drell- 
Yan, 11^-boson and jet production in hadronic collisions at Fermilab, sophisticated global fitting 



efforts were developed [4, 5, 6, 7] tfiat now include perturbative corrections calculated to next- 
to- leading order (or higher) in the strong coupling constant ag. With the increasing energies 
available at facilities such as the DESY ep collider HERA (and in future the Large Hadron 
Collider at CERN), PDF studies turned their attention primarily to exploring the region of very 
small X (down to x ~ 10~^), where the structure of the nucleon is dominated by its sea quark 
and gluon distributions. However, while opening the door to exploration of phenomena such 
as saturation and evolution in new kinematic regimes, one can argue whether DIS at low x 
measures the intrinsic structure of the nucleon or the hadronic structure of the virtual photon, 
7*. Because the virtual photon in the DIS process can fluctuate into qq pairs whose coherence 
length A ~ 1/Mx becomes large at small x, DIS at x<0.1 really probes the 7* A'' interaction 
rather than the structure of the nucleon or the 7* separately. At high x, in contrast, the virtual 
photon is point-like and unambiguously probes the structure of the nucleon [8] . 

Moreover, despite the impressive achievements over the past 4 decades, there are still 
some regions of kinematics where our knowledge of structure functions and PDFs remains 
unacceptably poor, with little progress made since the 1970's. A striking example is the region 
of large x (x>0.5), where most of the momentum is carried by valence quarks, and sea quarks 
and gluons are suppressed. Here the valence quark PDFs can be more directly related to quark 
models of hadron structure; however, the rapidly falling cross sections have made precision 
measurements extremely challenging. Another example is the pre-asymptotic region dominated 
by nucleon resonances, where data on the individual transverse and longitudinal cross sections at 
intermediate and high values of are either nonexistent or have large uncertainties. With the 
availability of continuous, high luminosity electron beams at the CEBAF accelerator, the first 
decade of experiments at Jefferson Lab has seen a wealth of high-quality data on unpolarized 
structure functions of the nucleon, penetrating into the relatively unexplored large-x domain 
and the transition region between resonances and scaling. 

The new data in the resonance region confirmed in spectacular fashion the phenomenon 
of quark-hadron duality in the proton F2 structure function, and revealed intriguing details 
about the workings of duality in a number of other observables. The impact of this has been 
a re-evaluation of the applicability of perturbative QCD to structure functions at low Q^, and 
has allowed a much larger body of data to be used in global PDF analyses [9]. Jefferson Lab 
has also set a new standard in the determination of Rosenbluth-separated longitudinal and 
transverse structure functions, which eliminates the need for model-dependent assumptions that 
have plagued previous extractions of structure functions from cross section data. 

On the theoretical front, the region of large x and low brings to the fore a number of issues 
which complicate structure function analysis, such as suppressed target mass and higher 

twist corrections, and nuclear corrections when scattering from nuclear targets. Controlling these 
corrections requires more sophisticated theoretical tools to be developed, and has motivated 
theoretical studies, many of which are still ongoing. It has also paved the way towards the 
12 GeV experimental program, in which structure functions will be measured to very high x in 
the DIS region, addressing some long-standing questions about the behavior of PDFs as x — )■ 1. 

In the next section we summarize the kinematics and formalism relevant for inclusive 
lepton-nucleon scattering, including the key results from the operator product expansion. 
Measurements of the proton structure functions and their moments are reviewed in Sec. 3, 
together with their role in the verification of quark-hadron duality. Data on the deuteron 
F2 structure function are presented in Sec. 4, and the extraction from these of the neutron 
F2 structure function F2 is discussed in Sec. 5. Section 6 reviews new measurements of 
the longitudinal structure function Fl, while Sec. 7 surveys results from semi-inclusive pion 
production. Finally, in Sec. 8 we describe the impact that the Jefferson Lab data have had 
on our understanding of nucleon structure in a global context, and briefly outline prospects for 
future measurements in the 12 GeV era. 



2. Formalism 

2.1. KinemMics 

Because of the small value of the electromagnetic fine structure constant, a = e^/Air, the 
inclusive scattering of an electron from a nucleon, e{k) + N{p) — >■ e'{k') + X, can usually be 
approximated by the exchange of a single virtual photon, 7* ((7), where q = k' — k. In terms 
of the laboratory frame incident electron energy E, the scattered electron energy E' , and the 
scattering angle 6, the photon virtuality is given hy -q^ = Q"^ = AEE'sni^e/2, where the 
electron mass has been neglected. The invariant mass squared of the final hadronic state X is 
VF^ = {p + q)"^ = M2 + 2Mz/-Q2 = M"^ + Q'^{l-x)/x, where M is the nucleon mass, v = E-E' 
is the energy transfer, and x = Q"^ /2Mu is the Bjorken scaling variable. 

In the one photon exchange approximation the spin-averaged cross section for inclusive 
electron-nucleon scattering in the laboratory frame can be written as 



dSldE' Q4 E 

where the leptonic tensor L^^ = 2 (jt^kl + k'Jtv — g^^k ■ k'^ , and using constraints from Lorentz, 
gauge and parity invariance the hadronic tensor W^'^ can in general be written as 

MH^- = + ^) Q^) + (p- - ^,^) (p^ - ^,^) . (3) 

The structure functions Fi and F2 are generally functions of two variables, but become 
independent of the scale in the Bjorken limit, in which both and v ^ 00 with x fixed. 
At finite values of a modified scaling variable is more appropriate [10, 11], 

e=-^, withp=M = ^l + Q2/^2 ^ (4) 

which tends to x in the Bjorken limit. 

In terms of cross sections for absorbing helicity ±1 (transverse) and helicity (longitudinal) 
photons, or and a^, the cross section can be written as 



(Pa 
dndE' 



r(aT{x,Q^) + eaL{x,Q^)), (5) 



where F = {a/2TT'^Q'^){E' /E)K/{1 — e) is the flux of transverse virtual photons, with the factor 
K = u{l — x) in the Hand convention [12], and 



1 + 2 




(6) 



is the relative flux of longitudinal virtual photons. Equating Eqs. (2) and (5), the structure 
functions can be written in terms of the photoabsorption cross sections as 

Fi(x,Q2) = -^M(7t(x,Q2), (7) 
= 4^^(i + j2/Q2) [rT{.x,Q^) + aL{x,Q^)\, (8) 

which reveals that F\ is related only to the transverse virtual photon coupling, while F2 is 
a combination of both transverse and longitudinal couplings. One can also define a purely 
longitudinal structure function Fx,, 



Fl = p^F2 - 2xFi = 2xFi R , 



(9) 



where R = aL/cfr is the ratio of longitudinal to transverse cross sections. 

The separation of the unpolarizcd structure functions into longitudinal and transverse parts 
from cross section measurements can be accomplished via the Rosenbluth, or longitudinal- 
transverse (LT), separation technique [13], by making measurements at two or more values of e 
for fixed x and Q^. Fitting the reduced cross section a/T linearly in e yields (and therefore 
Fi) as the intercept, while the ratio R is obtained from the slope. Note that F2 can only be 
extracted from cross sections either by measuring at e = 1 or by performing LT separations. At 
typical Jefferson Lab kinematics the contribution of Fl to F2 can be significant. 

The above discussion assumed the dominance of the one-photon exchange amplitude in 
describing the neutron current electron scattering cross section. In principle there are additional 
contributions arising from the exchange of a Z boson, and in particular the interference 
between 7* and Z exchange. The interference is in fact very relevant for parity-violating 
electron scattering, discussed elsewhere in this volume in connection with extractions of strange 
electromagnetic form factors from parity- violating asymmetries. 



2.2. Operator product expansion 

The theoretical basis for describing the dependence of structure functions in QCD is Wilson's 
operator product expansion (OPE) [14]. The quantities most directly amenable to a QCD 
analysis are the moments of structure functions, the n-th moments of which are defined as 

mJ")(Q2) = /'dxx"-^Fi(x,Q2), M^%Q^) = [\xx''-^F2,l{x,Q^). (10) 
Jo ' Jo 

As will become relevant in the discussion of duality in Sec. 3 below, at large ^ ^qcd 
the moments can be expanded in powers of with the coefficients in the expansion given 

by matrix elements of local operators corresponding to a certain twist, r, defined as the mass 
dimension minus the spin, n, of the operator. For the n-th moment of F2, for instance, one has 
the expansion 

^(n)^^2) ^ g „ = 2,4,6... (11) 

T=2,4... ^ 

where ^t"^ are the matrix elements with twist < r. As the argument suggests, the dependence 

(n) 

of the matrix elements can be calculated perturbatively, with A)- expressed as a power series 
in asiQ"^). For twist two, the coefficients are given in terms of matrix elements of spin-n 

operators, Ag"^'*^ ' ' H = {p\'ipj^'^^ iDi^^ ■ ■ ■ iDi^"^ tplp) , where is the quark field, is 

the covariant derivative, and the braces {• • •} denote symmetrization of indices and subtraction 
of traces. 

The leading-twist terms correspond to diagrams such as in Fig. 1 (a), in which the virtual 
photon scatters incoherently from a single parton. The higher-twist terms in Eq. (11) are 
proportional to higher powers of l/Q^ whose coefficients are matrix elements of local operators 
involving multi-quark or quark-gluon fields, such as those depicted in Fig. 1 (b) and (c). The 
higher twists therefore parametrize long-distance multi-parton correlations, which can provide 
clues to the dynamics of quark confinement. 

The additional terms (referred to as the "trace terms") in the twist- two matrix elements 
involve structures such as p^^'^'W and are thus suppressed by powers of p^/Q^ ~ Q'^/v^. While 
negligible in the Bjorken limit, at finite these give rise to the so-called target mass corrections 
(TMCs), and arc important in the analysis of Jefferson Lab data at large values of x. Because 
their origin is in the same twist-two operators that give rise to structure function scaling, TMCs 




Figure 1. (a) Leading- twist ("handbag") contribution to the structure function, (b) higher- 
twist ("cat's ears") four-quark contributions, (c) higher-twist quark-gluon interactions. 



are formally twist-two effects and are kinematical in origin. Inverting the expressions for the full 
moments including the trace terms, the resulting target mass corrected F2 structure function in 
the OPE is given by [15, 16] 



pTMC 



(x,Q2 



2^3-' 2 



+ 



du 1 



2M'^x 



(12) 



where F^^^ is the structure function in the M"^ jQ^ limit. Similar expressions are found for 

the F\ and F]^ structure functions [15, 16]. One should note, however, that the OPE result for 
TMCs to structure functions is not unique; in the collinear factorization approach, for example, 
in which parton distributions are formulated a "priori in momentum space, different expressions 
for TMCs arise [17]. While both formalisms give the same results in the Bjorken limit, the 
differences between these at finite can be seen as representing an inherent prescription 
dependence and systematic uncertainty in the analysis of structure functions at low . 



3. Proton Fi structure function 

Measurements of the proton F^ structure function have been taken at Jefferson Lab over a 
range of kinematics, from as low as 0.1 GeV^ and below (to study the transition to the 
photoproduction point, = 0) and up to = 8 GeV^ (to study the large-x behavior 
and quark-hadron duality). At the larger values the high luminosity provided by the 
CEBAF accelerator has allowed significant improvement in the statistical precision of high- 
X measurements over all previous experiments. In addition, with the HMS spectrometer in 
Hall C well understood, LT separated cross sections have been measured with better than 1.6% 
systematic point-to-point uncertainties and typically less than 1.8% normalization uncertainties. 

Precision F^ spectra extracted from cross sections measured in Hall C [18, 19, 20] are shown 
in Fig. 2 (left) as a function of x for several values (Q^ = 0.5, 1.5, 3 and 5.5 GeV^), together 
with previous SLAG [21] and NMC [22] data at lower x. The data have been bin-centered to 
the common values shown for all measurements within the range of 20% of the central value, 
utilizing a fit [23] to the DIS data and a global fit [24] to Jefferson Lab resonance region data. 
In addition to the Hall G measurements, there now also exists a large body of F^ data from 
Hall B covering a significant range of kinematics, which is afforded by the large acceptance of 
the GLAS spectrometer. An example of the F2 spectrum extracted from GLAS is shown in 
Fig. 2 (right) at = 0.775 GeV^. In Table 1 we present a complete list of all unpolarized 
inclusive and semi-inclusive measurements on the proton and deuteron performed at Jefferson 
Lab through 2010, including their current status. 



3.1. Quark-hadron duality 

The proton F2 data in Fig. 2 illustrate the intriguing phenomenon of quark-hadron duality, 
which relates structure functions in the nucleon resonance and DIS regions. First observed by 




Figure 2. (Left) Proton structure function data from Hall C [18, 19, 20], SLAC [21], 
and NMC [22] at = 0.5, 1.5,3 and 5.5 GeV^, compared with a fit [24] to the transverse and 
longitudinal resonance cross sections from photoproduction to = 9 GeV^ (solid), and a global 
fit [23] to DIS data (dashed). (Right) Proton Ff from CLAS in Hall B at = o.775 GeV^ 
(stars) [25] , compared to earlier Hall C data (open circles) [26] . 



Bloom and Oilman [37] in the early inclusive SLAC data (hence also referred to as "Bloom- 
Gilman duality" ) , the structure functions in the resonance region are found on average to equal 
the structure functions measured in the "scaling" region at higher W. The resonance data 
oscillate around the scaling curve and slide along it with increasing Q^, as seen in Fig. 2 (left). 

The early Fg data from SLAC were extracted from cross sections assuming a fixed value 
for R (= 0.18), and with a scaling curve parametrizing the limited data available in the early 
1970's [38]. Since the original measurements, the Fg structure function has become one of 
the best studied quantities in lepton scattering, with data from laboratories around the world 
contributing to a global data base spanning over five orders of magnitude in both x and . With 
the advent of the Jefferson Lab data, precise F2 measurements now also exist in the resonance 
region up 

to Q2 ~ 8 GeV^, allowing many new aspects of duality to be quantified for the first 

time [39]. 

While the early duality studies considered only the F2 structure function [37], Jefferson Lab 
experiments have in addition revealed the presence of duality in other observables. For example. 
Fig. 3 shows new LT-separated data from Jefferson Lab experiment E94-110 for the proton 
transverse (Ff) and longitudinal (F^) structure functions in the nucleon resonance region [18]. 
LT-separated data from SLAC, predominantly in the DIS region, are also shown for comparison 
[42]. Where they refer to the same kinematic values, the Jefferson Lab and SLAC data are 
in excellent agreement, providing confidence in the achievement of the demanding precision 
required of this type of experiment. In all cases, the resonance and DIS data merge smoothly 
with one another in both x and Q^. 

The availability of leading twist PDF-based global fits [4, 5, 6, 7, 9] allows comparison of 
the resonance region data with leading twist structure functions at the same x and Q^. The 
resonance data on the Fi and Fl structure functions are also found to oscillate around the 
perturbative QCD (pQCD) curves, down to as low as 0.7 GeV^. Because most of the 
data lie at large values of x and small Q^, it is vital for tests of duality to account for the 
effects of kinematical target mass corrections [15], which give large contributions as x — )> 1 [16]. 



Table 1. Listing of Jefferson Lab unpolarized inclusive and semi-inclusive electron- nucleon 
scattering experiments. 



Experiment Hall Target Observable 



Reference Status 



E94-110 

E99-118 

CLAS 
EOO-002 

EOO-108 
EOO-116 

E02-109 

E03-012 
(BoNuS) 

E06-009 



C 



B 
C 

C 
C 

C 
B 

C 



p R in resonance region 

p, d nuclear dependence of 
R at low 

d inclusive cross sections 

p, d F2 at low 



p, d semi-inclusive tt^ 
electroproduction 

p, d inclusive resonance 
region cross sections 
at intermediate 

d R in resonance region 



d{n) neutron F2 via 
spectator tagging 



R in resonance region 

& beyond: extension of 
E02-109 to Q2 = 4 GeV^ 



[18, 27, 28, 29] data taken in 1999, 
analysis completed 

[30] data taken in 2000, 

analysis completed 

[25, 31] el/e2 run periods 

[20] data taken in 2003, 

analysis completed, 
pub. in progress 

[32, 33] data taken in 2003, 
analysis completed 

[19] data taken in 2003, 

analysis completed 



[34] data taken in 2005, 

analysis in progress 

[35] data taken in 2005, 

analysis completed, 
pub. in progress 

[36] data taken in 2007, 

analysis in progress 



This is clear from Fig. 3, where the data are compared with leading twist structure functions 
computed from PDFs to ncxt-to-ncxt-to-lcading order (NNLO) accuracy from Alckhin [40] and 
MRST [41, 43]. The latter are shown with (solid) and without (dotted) target mass corrections, 
and clearly demonstrate the importance of subleading effects at large x. In particular, 

TMCs give additional strength at large x observed in the data, which would be significantly 
underestimated by the leading twist functions without TMCs. 

The phenomenological results raise the question of how can a scaling structure function 
be built up entirely from resonances, each of whose contribution falls rapidly with [44]? A 
number of studies using various models have demonstrated how sums over resonances can indeed 
yield a independent function (see Ref. [39] for a review). The key observation is that while 



X 



X 



Figure 3. Purely transverse (left) and longitudinal (right) proton structure functions 2xFi and 
Fl in tiie resonance region [18] (triangles), compared with earlier data from SLAC (squares). 
The curves are leading twist structure functions computed at NNLO from Alekhin (dashed) 
[40] and MRST [41] with (solid) and without (dotted) target mass corrections. The prominent 
resonance regions (A, ^n, F15) are indicated by the arrows along the abscissa. 



the contribution from each individual resonance diminishes with Q^, with increasing energy 
new states become accessible whose contributions compensate in such a way as to maintain 
an approximately constant strength overall. At a more microscopic level, the critical aspect 
of realizing the suppression of the higher twists is that at least one complete set of even 
and odd parity resonances must be summed over for duality to hold [45]. For an explicit 
demonstration of how this cancellation takes place in the SU(6) quark model and its extensions, 
see Refs. [45, 46, 47]. 

3.2. Structure function moments 

The degree to which quark-hadron duality holds can be more precisely quantified by computing 
integrals of the structure functions over x in the resonance region at fixed values, 
^xtC '^^ P^ix-jQ"^), where Xth corresponds to the pion production threshold at fixed Q^, and 
Xres = Q^/ (W^res~-^^+Q^) indicates the x value at the same where the traditional delineation 
between the resonance and DIS regions aiW = Wres = 2 GeV is made. These integrals can then 
be compared to the corresponding integrals of the structure functions fitted to the higher-VF, 
deep-inelastic data, at the same and over the same interval of x. The early phenomenological 
findings [26] suggested that the integrated strength of the resonance structure functions above 

~ 1 GeV^ was indeed very similar to that in the deep-inelastic region, including in each of 
the individual prominent resonance regions. In this section we explore the duality between the 
resonance and deep-inelastic structure functions in the context of QCD moments. 

According to De Rujula, Georgi and Politzer [48], one can formally relate the appearance of 
quark-hadron duality to the vanishing suppression of higher twist matrix elements in the QCD 
moments of the structure functions [14]. Namely, if certain moments of structure functions 
are observed to be independent of Q^, as implied by duality, then from Eq. (11) the moments 
must be dominated by the leading, independent term, with the 1/Q'^~'^ higher twist terms 
suppressed. Duality is then synonymous with the suppression of higher twists, which in partonic 
language corresponds to the suppression, or cancellation, of interactions between the scattered 



quark and the spectator system such as those illustrated in Fig. 1 (b) and (c). Conversely, if 
the moments display power-law dependence, then this implies violation of duality; moreover, 
if the violation is not overwhelming, the dependence of the data can be used to extract 
information on the higher twist matrix elements [49]. 
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Figure 4. Total n = 2 moments of the proton i^f (top), Ff (center) and Ff^ (bottom) structure 
functions determined from global fits to existing DIS data and Jefferson Lab resonance region 
data [50] , compared with moments computed from the leading twist PDFs from MRST at NNLO 
[41]. 



The first determination [51] of the F2 moments from Jefferson Lab data was made utilizing 
structure functions measured in Hall C [26], while a later evaluation [25] included the large body 
of additional data from CLAS. More recently, the extraction of F2 has been further enhanced 
by LT-separated data from Hall C, shown in Fig. 2 along with the fit [24] to the LT separated 
cross sections. The n = 2 moments for the proton Fg, -Ff , and F£ structure functions are shown 
in Fig. 4 versus Q^, as determined from integrating this fit. For they are found to be in very 
good agreement with the earlier measurements. Also shown is the leading-twist contribution 
calculated from the MRST parameterization [41], corrected for target mass effects [15]. 

One of the most striking features of the results in Fig. 4 is that the elastic-subtracted moments 
exhibit the same dependence as the PDF fits down to ~ 1 GeV^. Even with the elastic 
contribution included, which vanishes in the Bjorken limit and is hence pure higher twist, there 
is excellent agreement between the resonance and DIS data for > 2 GeV^. Until very recently 
[6, 9], this fact has not been widely appreciated or utilized in global PDF fitting efforts [4, 5], 
which typically impose cuts on data of > 4 GeV^ and W"^ > 12 GeV^. 

While the OPE provides a systematic approach to identifying and classifying higher twists, 
it does not reveal why these are small or how duality is realized globally and locally. To further 
explore the local aspects of duality within a perturbative QCD context, a ground-breaking new 



approach using "truncated" moments of structure functions, developed recently by Forte et al. 
[52] and extended by Kotlorz & Kotlorz [53], was applied by Psaker et al. [54] to Jefferson Lab 
data. The virtue of truncated moments is that they obey a similar set of evolution equations 
as those for PDFs themselves, which therefore enables a rigorous connection to be made between 
local duality and QCD. It allows one to quantify the higher twist content of various resonance 
regions, and determine the degree to which individual resonances are dominated by the leading 
twist components. 

Defining the n-th truncated moment Ain of a PDF q{x, Q^) between Xmin and Xmax as 



The symbol (8> denotes the Mellin convolution of the truncated moment and the "splitting 
function" P^, which is related to the usual DGLAP evolution splitting function P [55] by 
Pli{z,as{Q'^)) = P{z,as{Q^)). The extent to which nucleon structure function data in 
specific regions in x (or W) are dominated by leading twist can be determined by constructing 
empirical truncated moments and evolving them to different Q^. Deviations of the evolved 
moments from the experimental data at the new then reveal any higher twist contributions 
in the original data. 

A next-to-leading order (NLO) analysis [54] of data on the proton F^' structure function from 
Jefferson Lab covering a range in from 1 GeV^ to ~ 6 GeV^ revealed intriguing behavior 
of the higher twists for different nucleon resonance regions. Assuming that data beyond a 
large enough (taken to be = Qq = 25 GeV^ in Ref. [54]) are dominated by leading twist, 
the truncated moments were computed at Qq and evolved to lower Q^. Note that the truncated 
moments are computed over the range Wth < W < Wmaxj where the Wth = M + m.^ is the 
inelastic threshold. At = 1 GeV^ this corresponds to the integration range Xmin <x< xth, 
where a:;th = [l + 171^,(171,, + 2M) /Q^] ~ o.78. 

The ratio of the truncated moments of the data to the leading twist is shown in Fig. 5(a) as 
a function of Wmax at = 1 GeV^, with and without target mass corrections. Including the 
effects of TMCs, the leading twist moment differs from the data by ~ 15% for Wmax > 1-5 GeV. 
To quantify the higher twist content of the specific resonance regions, and at different values of 

, several intervals in W are considered: W'^^^ < < 1.9 GeV^ (A(1232) or first resonance 
region); 1.9 <W'^ < 2.5 GeV^ (S'ii(1535) or second resonance region); and 2.5 <W'^ < 3.1 GeV^ 
(Fi5(1680) or third resonance region). The higher twist contributions to in these regions 
are shown in Fig. 5(b) as ratios to moments of the data. 

The results indicate deviations from leading twist behavior of the entire resonance region data 
(filled circles in Fig. 5(b)) at the level of <15% for all values of considered, with significant 
dependence for Q^<4 GeV^. The strong dependence of the higher twists is evident 
here in the change of sign around = 2 GcV^, with the higher twists going from ~ —10% at 
Q2 = 1 GeV^ to ^ 10-15% for Q2 ^ 3 _ 6 GeV^. At larger the higher twists are naturally 
expected to decrease, once the leading twist component of the moments begins to dominate. 

Interestingly, the magnitude of the higher twist contributions in the A region (diamonds) is 
smallest, decreasing from « —15% of the data at = 1 GeV^ to values consistent with zero 
at larger Q^. The higher twists are largest, on the other hand, for the Su region (squares), 
where they vary between « -15% of the data at = 1 QeV^ and 20-25% at ~ 5 GeV^. 
Combined, the higher twist contribution from the first two resonance regions (dotted curve) is 




(13) 



the evolution equations for the truncated moments can be written as 




(14) 




Figure 5. (a) Ratio of the M.2 truncated moments of the data to the leading twist + TMC 

(sohd), and data to leading twist without TMC (dashed) at = 1 GcV^, as a function of Vl^max- 
(b) dependence of the fractional higher twist (HT) contribution to the n = 2 truncated 
moment data, for various intervals in W. 



< 15% in magnitude for all Q^. The rather dramatic difference between the A and the Su, may, 
at least in part, be due to the choice of the differentiation point of W"^ = 1.9 GeV^. A lower 

choice, for instance, would lower the higher twist content of the at large Q^, while raising 
that of the A. However, this W'^ choice corresponds to the local minimum between these two 
resonances in the inclusive spectra, and is the one most widely Titilized. 

The higher twist content of the Fi^ region (open circles) is similar to the Su at low Q^, but 
decreases more rapidly for > 3 GeV^. The higher twist content of the first three resonance 
regions combined (dashed curve) is < 15-20% in magnitude for < 6 GeV^. Integrating up 
to H^max = 4 GeV^ (filled circles), the data on the n = 2 truncated moment are found to be 
leading twist dominated at the level of 85-90% over the entire range. 

The results in Fig. 5 can be compared with quark model expectations [45, 46], which predict 
systematic deviations of resonance data from local duality. Assuming dominance of magnetic 
coupling, the proton data are expected to overestimate the DIS function in the second and 
third resonance regions due to the relative strengths of couplings to odd-parity resonances; the 
positive higher twists observed in Fig. 5(b) for Q^>2 GeV^ indeed support these predictions. 

4. Deuteron F2 measurements 

Together with hydrogen, inclusive lepton scattering from deuterium targets has provided an 
extensive data base of F2 structure function measurements over a large range of kinematics. 
While a significant quantity of F2 data was collected from experiments at CERN and SLAG, 
the quasi-elastic and nucleon resonance regions, especially at low and moderate (~ few 
GeV^), were only mapped precisely with the advent of Jefferson Lab data [26, 31]. Similarly, 
Jefferson Lab contributed with precision F2 data in the region of x > 1 [56], of relevance for 
constructing higher moments of deuteron structure functions. 

Experiments in Hall C [26, 56] have provided F2 data in select regions of x and Q"^, while 
inclusive cross section measurements in GLAS have covered a continuous two-dimensional region 
over the entire resonance region up to = 6 GeV^. This combination is rather useful for 
determining moments of F2. In such extractions one usually assumes that the ratio R for the 
deuteron is similar to that for the proton at scales of a few GeV^ [57]. New measurements 



which will test this assumption will be reviewed in Sec. 6. 

Table 2. Lowest two moments (for n = 2 and 4) of the isovector F2 structure function. 
Experimental results for 4 GeV^ are compared with lattice calculations extrapolated to 

the chiral limit. 



n 


Niculescu et al. [58] 


Osipenko et al. [59] 


Detmold et al. [60] 




(Hall C) 


(Hah B) 


(lattice) 


2 


0.049(17) 


0.050(9) 


0.059(8) 


4 


0.015(3) 


0.0094(16) 


0.008(3) 



The results of the moment analyses are shown in Table 2, expressed as the isovector (proton 
minus neutron) combination, and compared with isovector moments from lattice QCD [60, 61]. 
For simplicity the neutron moments here are defined as the difference between the deuteron and 
proton moments — see, however. Sec. 5.1 below. The n = 2 moments from the Hall B [59] and 
Hall C [58] analyses agree well with each other, and with the lattice extraction, which includes 
the effects of pion loops and the intermediate A(1232) resonance in the chiral extrapolation. 
For the n = 4 moment the comparison between the Hall B and Halls C results shows a slight 
discrepancy, which may be reduced once precision high-x data at ~ 4 GeV^ are included in 
the extractions. 

Analysis of the Q^-dependcncc of the deuteron F2 moments in Ref. [31] suggests a partial 
cancellation of different higher twist contributions entering in the OPE with different signs, 
which is one of the manifestations of quark-hadron duality. The slow variation with of the 
structure function moments, down to ~ 1 GeV^, was also found in the analysis of proton 
data [51], where such cancellations were found to be mainly driven by the elastic contribution. 
Furthermore, by comparing the proton and (nuclear corrected) neutron structure function 
moments, the higher twist contributions were found to be essentially isospin independent [62]. 
This suggests the possible dominance of ud correlations over uu and dd in the nucleon, and 
implies higher twist corrections that are consistent with zero in the isovector F2 structure 
function. 

More recently, high precision deuteron cross sections in the resonance region have been 
measured in Hall C [34, 36] with the aim of providing LT separated deuteron structure functions 
of comparable precision and kinematic coverage to those performed for the proton. These higher 
precision LT separated data will allow for a significant further reduction in the uncertainties in 
the current nonsinglet F2 extractions. 

5. Neutron F2 structure function 

A complete understanding of the valence quark structure of the nucleon requires knowledge of 

both its u and d quark distributions. While the u distribution is relatively well constrained 
by measurements of the proton F2 structure function, in contrast the d quark distribution is 
poorly determined due to the lack of comparable data on the neutron structure function F^*. 
The absence of free neutron targets makes it necessary to use light nuclei such as deuterium as 
effective neutron targets, and one must therefore deal with the problem of extracting neutron 
information from nuclear data. 

5.1. Neutron structure from inclusive F2 data 

In standard global PDF analyses, sensitivity to the d-quark from charged lepton scattering is 
primarily provided by the neutron in the deuteron. Usually the neutron F2 structure function 



is extracted by subtracting the deuteron and proton structure function data assuming that 
nuclear corrections are neghgible. At large x, however, the ratio of the deuteron to free nucleon 
structure functions is predicted to deviate significantly from unity [63, 64, 65, 66], which can 
have significant impact on the behavior of the extracted neutron structure function at large x 
[9, 67]. 

Even when nuclear effects are considered, there exist practical difficulties with extracting 
information on the free neutron from nuclear data, especially in the nucleon resonance region, 
where resonance structure is largely smeared out by nucleon Fermi motion. A recent analysis 
[68] used a new method [66] to extract F2 from F2 and F2 data, in which nuclear effects are 
parameterized via an additive correction to the free nucleon structure functions, in contrast to 
the more common multiplicative method [69] which fails for functions with zeros or with non- 
smooth data. In the standard impulse approximation approach to nuclear structure functions, 
the deuteron structure function can be written as a convolution [65, 66] 

Fiix,Q^) = J2 ldyfN/d{y,p)Fi'ix/y,Q^), (15) 

N=p,n 

where fN/d is the light-cone momentum distribution of nucleons in the deuteron (or "smearing 
function" ) , and is a function of the momentum fraction y of the deuteron carried by the struck 
nucleon, and of the virtual photon "velocity" p (see Eq. (4)). The smearing function encodes the 
effects of the deuteron wave function, accounting for nuclear Fermi motion and binding effects, 
as well as kinematical finite-Q^ corrections. Although not well constrained, nucleon off-shell 
effects have also been studied [63, 64, 65]; their influence appears to be small compared with 
the errors on the existing data, except at very large x. 




Figure 6. (Left) Neutron F2 structure function extracted from inclusive deuteron and proton 
data at = 1.7 GeV^ [68], together with the reconstructed F2. (Right) Neutron structure 
function extracted from the BoNuS experiment [35], compared with the Bosted/Christy model 
of the neutron F2 [70] and the corresponding Christy /Bosted parametrization for F2 [24]. 

In Fig. 6 we illustrate the results of a typical extraction of F2 from Jefferson Lab F2 and 
F2 data at = 1.7 GeV^. The proton data show clear resonant structure at large x, which is 
mostly washed out in the deuteron data. The resulting neutron F2 is shown after two iterations 
of the procedure with an initial guess of F2 = -Ff ■ Clear neutron resonance structure is visible 
in the first (A) and second resonance regions, at x 0.75 and 0.55, respectively, with some 



structure visible also in the third resonance region at x ~ 0.45, albeit with larger errors. The 
dcutcron F2 reconstructed from the proton and extracted neutron data via Eq. (15) indicates 
the relative accuracy and self-consistency of the extracted ¥2- 

Of course it is not possible to avoid nuclear model dependence in the inversion procedure, 
and some differences in the extracted F2 will arise using different models for the smearing 
functions /jv/d- To remove, or at least minimize, the model dependence in the extracted free 
neutron structure, several methods have been proposed, such as utilizing inclusive DIS from 
v4 = 3 mirror nuclei [71, 72, 73], and semi-inclusive DIS from a deuteron with spectator tagging 
[35, 74]. In addition, experiments involving weak interaction probes [75, 76, 77] can provide the 
information on the flavor separated valence quark distributions directly. In the next section we 
discuss in detail one of these methods for determining the free neutron structure, namely the 
BoNuS experiment at Jefferson Lab [35]. 

5.2. Tagged neutron structure functions 

To overcome the absence of free neutron targets, the Hall B BoNuS (Barely Off-shell NeUtron 
Structure) experiment has measured inclusive electron scattering on an almost free neutron using 
the CLAS spectrometer and a recoil detector to tag low momentum protons. The protons are 
tracked in a novel radial time projection chamber [78] utilizing gas electron multiplier foils to 
amplify the proton ionization in a cylindrical drift region filled with a mixture of Helium and 
Di-Methyl Ether, and with the proton momentum determined from the track curvature in a 
solenoidal magnetic field. Slow backward-moving spectator protons are tagged with momenta 
as low as 70 MeV in coincidence with the scattered electron in the reaction D{e,e'ps)X. This 
technique ensures that the electron scattering took place on an almost free neutron, with its 
initial four-momentum inferred from the observed spectator proton. Cutting on spectator 
protons with momenta between 70 and 120 MeV and laboratory angles greater than 120 degrees 
minimizes the contributions from final state interactions and off-shell effects to less than several 
percent on the extracted neutron cross section. 

While inclusive scattering from deuterium results in resonances which are significantly 
broadened (often to the point of being unobservable) , determination of the initial neutron 
momentum allows for a dramatic reduction in the fermi smearing effects and results in 
reconstructed resonance widths comparable to the inclusive proton measurements. In addition, 
the large CLAS acceptance for the scattered electron allowed for the tagged neutron cross section 
to be measured over a significant kinematic range in both W'^ and at beam energies of 4.2 GeV 
and 5.3 GeV. BoNuS F2 data extracted at a beam energy of 5.3 GeV and = 1.7 GeV^ are 
shown in Fig. 6 (right panel) for x values from pion production threshold through the resonance 
region and into the DIS regime. This will allow for the first time a unambiguous study of the 
inclusive neutron resonance structure functions. 

An extension of BoNuS has been approved [79] to run at a beam energy of 11 GeV after the 
energy upgrade of the CEBAF accelerator. The kinematic coverage will allow the extraction 
of the ratio of neutron to proton structure functions F2/F2 to x as large as 0.8, and the 
corresponding d to u large x parton distribution ratio. Other quantities which BoNuS will 
make possible to measure include the elastic neutron form factor, quark-hadron duality on the 
neutron, semi-inclusive DIS and resonance production channels, hard exclusive reactions such 
as deeply- virtual Compton scattering or deeply- virtual meson-production from the neutron, as 
well as potentially the inclusive structure function of a virtual pion. 

6. Longitudinal structure function Fl 

The unpolarized inclusive proton cross section contains two independent structure functions. 

While F2 has been measured to high precision over many orders of magnitude in both x 
and Q^, measurements of the longitudinal structure function F^ (and the ratio R) have been 



significantly more limited in both precision and kinematic coverage. This is due in part to 

the challenges inherent in performing LT separations, which typically require point-to-point 
systematic uncertainties in e to be smaller than 2% to obtain uncertainties on Fl of less than 
20%. 

While the inclusive cross section is proportional to 2xFi + (F^, the F2 structure function 
~ 2xFi + F]^. so that F2 is onlv proportional to the cross section for e = 1. At high the 
scattering of longitudinal photons from spin- 1/2 quarks is suppressed, and in the parton model 
one expects Fl (and R) to vanish as — )• 00. At low Q'^, however, Fl is no longer suppressed, 
and could be sizable, especially in the resonance region and at large x. On the other hand, Fl 
is dominated by the gluon contribution at small x, where new measurements from HERA [80] 
have shown that it continues to rise. In the kinematic range of Jefferson Lab Fl has been found 
to be typically 20% of the magnitude of F2, which is consistent with earlier SLAC measurements 
where the kinematic regions overlap. 

An extensive program of LT separations has been carried out in Hall C, including 
measurements of the longitudinal strength in the resonance region for 0.3 < < 4.5 GeV^ 
for both proton [18] and deuteron [34, 36] targets. The Jefferson Lab experiments listed in 
Table 1 for which LT separations have been performed for the proton or deuteron are E94-110, 
E99-118, EOO-002, E02-109, and E06-009. The Jefferson Lab data complement well the previous 
results at smaller x from SLAC and NMC in this region, and improve dramatically on the 
few measurements that existed below = 8 GeV^ in this x region, which had typical errors on 
R and Fl of 100% or more — see Fig. 7 (left). 

The recent precision LT separated measurements of proton cross sections [18] have allowed 
for the first time detailed duality studies in all of the unpolarized structure functions and their 
moments. The results of the proton separated structure functions in the resonance region were 
presented in Fig. 3 in Sec.3. Although significant resonant strength is observed in Fl (or R), 
evidence of duality is nonetheless observed in this structure function, along with F2 and Fi . 

In addition to the proton data, Jefferson Lab experiment E02-109 measured the LT separated 
F2 and Fl structure functions of the deuteron, in the same and ranges, and with the 
same high precision as E94-110 did for the proton. This will allow quantitative studies of duality 
in both the longitudinal and transverse channels for the deuteron. If duality holds well for both 
the proton and neutron separately, it will hold to even better accuracy for the deuteron since the 
Fermi motion effects intrinsically perform some of the averaging over the resonances. However, if 
duality does not hold for the LT separated neutron structure functions, this should be observable 
in the deuteron data, and will thus provide a critical test for models of duality. 

In addition to the resonance region, measurements of the inclusive longitudinal proton and 
deuteron structure have also been performed at lower x (higher W"^) and lower Q^. While the 
longitudinal strength is significant at of several GeV^, the proton Fl structure function is 
constrained by current conservation to behave, for fixed W, as Fl ~ Q'^ for — > 0. However, 
even with the new Jefferson Lab data, which extend down to = 0.15 GeV^ [81], the at 
which this behavior sets in has not yet been observed. 

Another interesting test provided by the E99-118 data is whether the relative longitudinal 
contribution to the cross section embodied in R is different in the deuteron and proton at these 
low values. While the higher data from SLAC and NMC exhibit no significant difference 
in the deuteron and proton R, the Jefferson Lab results shown in Fig. 7 (right) suggest a 
possible suppression of R in the deuteron relative to the proton for < 1 GeV^. Although 
this suppression is consistent with the two lowest data points from SLAC, the uncertainties 
are dominated by systematic errors and the combined significance of the effect is still less than 
2 a. Conclusive experimental evidence for the possible suppression of R in deuterium at low 
will likely be provided when the results from additional data from EOO-002 are finalized in very 
near future. 




Figure 7. (Left) Sample LT separations from Jefferson Lab experiment E94-110 [18]. (Right) 
Difference between the ratios R in deuterium and hydrogen versus from E99-118 [81], 
compared with previous NMC and SLAC measurements. 



7. Semi-inclusive deep inelastic scattering 

In addition to the traditional F2 and Fl structure function observables, a number of other 
processes have been studied at Jefferson Lab over the past decade, with potentially important 
consequences for our understanding of the workings of QCD at low energy. In this section we 
focus on semi-inclusive pion electroproduction. 

An analysis of the semi-inclusive process eN — >■ ehX, where a hadron h is detected in 
the final state in coincidence with the scattered electron, has recently been made using data 
from Hall C in the resonance-scaling transition region [32, 33]. One of the main motivations 
for studying semi-inclusive meson production is the promise of flavor separation via tagging of 
specific mesons in the final state. In the valence quark region a produced vr"*" (tt") meson, for 
example, primarily results from scattering off a u (d) quark in the proton. 

The semi-inclusive cross section at LO in is given by a simple product of quark distribution 
and quark — > hadron fragmentation functions, 

^ ~ Y^elq{x)D^^{z) ^ M^{x,z). (16) 

Here the fragmentation function D^{z) gives the probability for a quark q to fragment to a hadron 
h with a fraction z = Ph ' p/q ■ P = E^/v of the quark's (or virtual photon's) laboratory frame 
energy. Although at LO the scattering and particle production mechanisms are independent, 
higher order pQCD corrections give rise to non-factorizable terms, which involve convolutions 
of the PDFs and fragmentation functions with hard coefficient functions [82] . 

For hadrons produced collinearly with the virtual photon, the invariant mass W of the 
undetected hadronic system X at large Q can be written [83] ^ + Q^{1- z){l - x)/x, 
where the hadron mass is neglected with respect to Q^. In the elastic limit, z — > 1, the hadron 
carries all of the photon's energy (with W' M), so z is also referred to as the "elasticity". 



While formally the LO factorized expression for the cross section (16) may be valid at 
large Q^, at finite there are important corrections arising from the finite masses of the 
target and produced hadron. One can show, however, that the LO factorization holds even 
at finite Q^, provided the parton distribution and fragmentation functions are expressed in 
terms of generalized scaling variables [84], q{x)Dl^{z) — )■ q{^h)Dq{Ch), where = {zh^/^x){l + 



^1 - Ax'^M'^mlJz'^Q'^) and Cft = C(l + fn\IC,uQ'^), with = ml + pl_^. Not surprisingly, 



these eff'ects become large at large x and z when is small; however, for heavier produced 
hadrons such as kaons or protons, significant effects can also arise at small values of z [84]. 

The validity of the factorized hypothesis in Eq. (16) relies on the existence of a sufficiently 
large gap in rapidity 77 = In [{Eh — pfi)/{Eh + p|)] /2 to allow a clean separation of the current 
fragmentation region (hadrons produced from the struck quark) from the target fragmentation 
region (hadrons produced from the spectator quark system). At high energies a gap of Arj « 2, 
which is typically required for a clean separation [85], can be achieved over a large range of z; at 
low energies, however, this can only be reached at larger values of z. On the other hand, at fixed 
X and the large-z region corresponds to resonance dominance of the undetected hadronic 
system X (corresponding to small W), so that the factorized description in terms of partonic 
distributions must eventually break down. It is vital therefore to establish empirically the limits 
beyond which the simple x and z factorization of Eq. (16) is no longer valid. 

It is intriguing in particular to observe whether W' can play a role analogous to W for duality 
in inclusive scattering, when the undetected hadronic system X is dominated by resonances 
W' < 2 GeV. In terms of hadronic variables the fragmentation process can be described through 
the excitation of nucleon resonances, A^*, and their subsequent decays into pions (or other 
mesons) and lower- lying resonances, N'*. The hadronic description must be rather elaborate, 
however, as the production of fast outgoing pions in the current fragmentation region at high 
energy requires nontrivial cancellations of the angular distributions from various decay channels 



where F^*Ar_j.Ar* is the N — >■ N* transition form factor, which depends on the masses of the 
virtual photon and excited nucleon (W = Mjv*), and T>N*^N'*h is a function representing the 

decay N* ^ N'*h. 

A dedicated experiment (EOO-018) to study duality in tt^ electroproduction was performed 
in Hall C [32, 33], in which a 5.5 GeV electron beam was scattered from proton and deuteron 
targets at between 1.8 and 6.0 GeV^, for 0.3 < x < 0.55, with z in the range 0.35 — 1. 
From the deuterium data the ratio of unfavored to favored fragmentation functions / 
was constructed, where corresponds to a pion containing the struck quark {e.g., 7r+ from 
a struck u or d quark), while D~ describes the fragmentation of a quark not contained in 
the valence structure of the pion {e.g., a d quark for the vr"*"). Since at moderate x the 
dependence on PDFs cancels, the fragmentation function ratio is approximately given by 



D-/D+ = (4 - Nl /N2 )/{AN2 /N^'^ - 1), where ATj is the yield of produced pions in 

Eq. (17). 



The Jefferson Lab data for D /D^ from EOO-108 are shown in Fig. 8 as a function of z at 
fixed X = 0.32 and = 2.3 GeV^ [32], and compared with earlier HERMES data at higher 
energies [87]. Despite the different energies, there is good overall agreement between the two 
measurements, even though the Jefferson Lab data sit slightly higher. Furthermore, the D~ /D~^ 
ratio extracted from the Jefferson Lab data shows a smooth dependence on z, which is quite 
remarkable given that the data cover the full resonance region, 0.88 < W'"^ < 4.2 GeV^. This 
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Figure 8. The ratio of unfavored to favored fragmentation functions D~ /D^ as a function of 
z extracted from deuterium data, for x = 0.32 [32]. 



strongly suggests a suppression or cancellation of the resonance excitations in the tt^/tt^ cross 
section ratio, and hence in the fragmentation function ratio. 

Similar cancellations between resonances naturally arises in quark models, such as those 
discussed by Close et al. [45, 47] for the tt^N'* reaction. The pattern of constructive and 

destructive interference, which was a crucial feature of the appearance of duality in inclusive 
structure functions, is also repeated in the semi-inclusive case when one sums over the states 
N'*. Moreover, the smooth behavior of the fragmentation function ratio D~ /D^ in Fig. 8 can 
be qualitatively understood from the relative weights of the matrix elements, which are always 4 
times larger than for vr" production. In this case the resonance contributions to this ratio cancel 
exactly, leaving behind only the smooth background as would be expected at high energies. This 
may account for the striking lack of resonance structure in the resonance region fragmentation 
functions in Fig. 8. 

8. Outlook 

8.1. Impact of Jefferson Lab data 

The first decade of unpolarized structure function measurements at Jefferson Lab has had 
significant impact on our understanding of nucleon structure, both for leading twist parton 
distributions and for the resonance-scaling transition and related studies of quark-hadron 
duality. With most of the data concentrated in the low-W region in the few GeV^ 

range, the greatest influence on the global data base has naturally been at large x. 

Recently a new global PDF analysis was performed [9], exploring the possibility of reducing 
the uncertainties at large x by relaxing the constraints on the kinematics over which data are 
included in the fit. The data sets combined proton and deuteron DIS structure functions from 
Jefferson Lab, SLAC and CERN (NMC) with new ep collider data from HERA, as well as new 
Drell-Yan, W asymmetry and jet cross sections from pp and pd collisions at Fermilab. The new fit 
(referred to as "CTEQ6X") allowed for a signiflcant increase in the large-x data set {e.g., a factor 
of two more DIS data points) by incorporating data for W'^ > 3 GeV^ and > 1.69 GeV^, 
lower than in the standard global fits [5, 88] which typically use W'^ > 12.25 GeV^ and 
> 4 GeV^. The new analysis also systematically studied the effects of target mass and 
higher twist contributions, and realistic nuclear corrections for deuterium data. 

Results from the CTEQ6X fit are shown in Fig. 9 (left) for the u and d quark PDFs, 
normalized to the earlier CTEQ6.1 fit, which had no nuclear or subleading corrections 
applied. The biggest change is the ~ 30-40% suppression of the d quark at x ~ 0.8, which is 
found to be stable with respect to variations in the W'^ and cuts, provided both TMC and 
higher twist corrections are included. The effect of the expanded data base is more dramatically 
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Figure 9. (Left) CTEQ6X fit for u and d quark PDFs, normalized to the earlier CTEQ6.1 fit 
[88]. The vertical lines show the approximate values of x above which PDFs are not directly 
constrained by data. The error bands correspond to Ax = 1. (Right) Relative errors on u and 
d quark PDFs, normalized to the relative errors in the reference fit. 

illustrated in Fig. 9 (right), which shows the relative u and d PDF errors for a range of W and 
cuts ("cutO" being the standard cut, "cut3" the more stringent CTEQ6X cut, and intermediate 
cuts "cutl" and "cut2"), normalized to those of a reference fit with "cutO" and no nuclear or 
subleading corrections. The result is a reduction of the errors by up to 40-60% at x>0.7, which 
will have a profound impact on applications of PDFs in high energy processes, such as those as 
the LHC, as well as in constraining low-energy models of quark distributions. 

8.2. Future prospects 

Uncertainties in PDFs will be further reduced with the availability of data at even larger x and 
from Jefferson Lab after its 12 GeV energy upgrade, which will determine the d quark distribution 
up to ~ 0.8 with minimal theoretical uncertainties associated with nuclear corrections. Several 
of the planned experiments include BoNuS12 [79], which will extend to larger x the earlier 
measurements of low-momentum, backward protons in semi-inclusive scattering from deuterium 
(Sec. 5.2); the MARATHON experiment [73], which plans to extract F2/F2 from the ratio of 
^He to structure functions, in which the nuclear corrections cancel to within ~ 1% [71, 72]; 
and the program of parity- violating DIS measurements on hydrogen [76] , which will be sensitive 
to a new combination of d/u in the proton, free of nuclear corrections. 

In the resonance region, experiment E12-10-002 [89] will extend proton and deuteron structure 
function measurements up to Q'^ ~ 17 GeV'^ and enable tests of quark-hadron duality over a 
much larger kinematic range, ultimately providing stronger constraints on large-x PDFs. And 
finally, a new avenue for exploring nucleon structure at 12 GeV will be opened up with semi- 
inclusive meson production experiments [90], which will test the factorization of scattering and 
fragmentation subprocesses needed for a partonic interpretation of semi- inclusive cross sections. 
A successful program of semi-inclusive measurements tagging specific mesons in the final state 
would allow unprecedented access to the flavor dependence of PDFs in previously unexplored 
regions of kinematics. 
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